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ABSTRACT 

There has been a recent interest in determining high statistical reliability in risk 
assessment of aircraft components. This report identifies the potential consequences of 
incorrectly assuming a particular statistical distribution for stress or strength data used in 
obtaining the high reliability values. The computation of the reliability is defined as the 
probability of the strength being greater than the stress over the range of stress values. 
This method is often referred to as the stress-strength model. 

A sensitivity analysis was performed involving a comparison of reliability results in 
order to evaluate the effects of assuming specific statistical distributions. Both known 
population distributions, and those that differed slightly from the known, were considered. 
Results showed substantial differences in reliability estimates even for almost nondetectable 
differences in the assumed distributions. These differences represent a potential problem 
in using the stress-strength model for high reliability computations, since in practice it is 
impossible to ever know the exact (population) distribution. 

An alternative reliability computation procedure is examined involving determination of 
a lower bound on the reliability values using extreme value distributions. This procedure 
reduces the possibility of obtaining nonconservative reliability estimates. Results indicated 
the method can provide conservative bounds when computing high reliability. 
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INTRODUCTION 


There has been an interest in quantitative reliability-based structural design for many 
years. An early example is the structural reliability development by Freudenthal. 1 Stress- 
strength reliability computations are a principal consideration in structural reliability design. 
Reliability methods have been considered for many structural applications including: civil 
engineering, 2 nuclear reactors, 3 fixed wing aircraft, 4 rotorcraft, 5 and space vehicle propulsion 
systems. 6 Very high structural reliability is expected to be achieved for most applications. A 
reliability goal of 0.9ra) per flight hour was suggested in 1955 by Lundberg 4 for fixed wing 
civil aircraft. Recently, Lincoln, 7,8 using reasoning similar to that of Lundberg, cited a reliabil- 
ity goal of 0.9(7) per flight for fixed wing military aircraft. The U.S. Army has instituted a 
new structural fatigue integrity criterion for rotorcraft which has been interpreted 5 as a re- 
quirement for a lifetime reliability of 0.9(6). 

The use of advanced materials whose structural properties are best characterized on a 
statistical basis appears to be a stimulant for increased interest in statistical-based structural 
design for airborne structures. 

A significant feature associated with predictions of structural reliability is that the 
consequence of a failure event may be more than reduced system performance or the 
inconvenience of a system being out of service; structural failure can be catastrophic in 
terms of loss of life and property. In this context it is imperative to evaluate the sensitiv- 
ity of structural reliability predictions to uncertainties. It appears that this issue has re- 
ceived little attention except for a brief note by Harris and Soms 9 and a recent 
presentation by Berens. 10 

There are many issues to be faced in obtaining quantitative structural reliability predic- 
tions. Such issues include system complexity (many components, multiple failure modes in 
each component, and interdependence of component behavior), sample or data set size 
associated with structural loading spectrum conditions and with mechanical properties, and 
the basis for characterizing structural qualification tests (the number of duplicate specimens 
and methods for compensation for untested effects such as the effect of environment). 

In addition, when predictions of structural behavior are required in the high reliability 
range, since sufficiently large data sets are usually not available, it is necessary to use 
parametric modeling methods. Assumed parametric functions permit extrapolation from 
available data to determine the probability of failure. Since the probability of failure is 
extremely small, this will always involve substantial extrapolation from what can be ob- 
served experimentally. The estimated reliability will therefore depend strongly on the as- 
sumed parametric probability density function (PDF). Slight deviation from the assumed 
model in tail regions can have a dramatic effect on high reliability estimates. 

In fact, one might argue, as does Freudenthal, 11 that because of the extrapolation in- 
volved, statistically-based high reliability calculations for complex systems must always be 
suspect: 

"When dealing with probabilities a clear distinction should be made between conditions 
arising in design of inexpensive mass products in which the probability figures are derived 
by statistical interpretation of actual observations or measurements (since a sufficiently 
large number of observations are actually obtainable), and conditions arising in design of 
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structures or complex systems. In the latter,, probability figures are used simply as a scale or 
measure of reliability that permits the comparison of alternative designs. The figures can 
never be checked by observations or measurement since they ale obtained by extrapolations 
so far beyond any possible range of observation that such extrapolation can no longer be 
based on statistical arguments but could only be justified by relevant physical reasoning. 

Under these conditions the absolute probability figures have no real significance — " 

Nonparametric stress-strength procedures do not require specific parametric assumptions, 
and so it mig ht be hoped that such procedures could circumvent this difficulty. However, 
Johnson 12 has noted that "The nonparametric approach has one serious drawback. In return 
for its distribution free property, it is not possible to establish high reliability even with mode- 
rate sample sizes." With respect to the use of parametric models, Box 3 has observed "all 
models are wrong, but some are useful," meaning that no parametric statistical model should 
be accepted uncritically. Whenever a model is used, it is the obligation of the analyst to in- 
vestigate the consequences of departures from an assumed model which, though small, are con- 
sistent with available data. Harris and Soms 9 has illustrated a "serious problem in the use of 
stress-strength relationships in estimating reliability." In particular, "stress-strength models in 
reliability theory are highly sensitive to small perturbations in extreme tails.” The perturba- 
tions considered may arise from an alternative mode of failure such as the presence of a flaw 
in a structure. Further, they note that the problem cannot be eliminated unless "astronomi- 
cally large sample sizes are employed.” 

In the following, the examination of the sensitivity of structural reliability estimates fo- 
cuses attention on one of the previously cited issues: the selection of a parametric PDF. 

The examination of the sensitivity of stress-strength reliability estimates is extended to addi- 
tional perturbation effects. The sensitivity of reliability estimates to the selection of paramet- 
ric models is considered with emphasis on graphical representations. The results are 
evaluated with regard to the usefulness of parametric stress-strength models for application to 
the high reliability regime of 0.9(6) to 0.9(7) when the consequence of failure may be cata- 
strophic. An alternative reliability computation procedure is examined involving determination 
of a lower bound on reliability which can be obtained independently of the assumed PDFs. 

STRESS-STRENGTH MODEL 

The statistical reliability as referred to in this report is determined in the following man- 
ner. Shown in Figure 1 is the stress-strength model where f 2 (s) and fi(S) represent the 
PDFs for the applied stress s and material strength S. 

Since the joint probability dR for the strength being greater than si can be written as, 

dR = f 2 (si) ds f°° fj(S)dS (1) 

*1 

then the reliability for all s values is 

R “-C few [£ f i( s ) ds ] d »- < 2 > 
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PROBABILITY DENSITY FUNCTIONS 

A wide variety of PDFs may be applied in obtaining R values. Some examples of PDFs 
are as follows: 

The PDF most often used in stress-strength models is the normal distribution (see Figure 1), 

W S ) = N( S ,^= 7 ^«p{-^(^) } , (3) 

where - » < S < », /i > 0, and a > 0. The mean of the population is /i, and the stan- 
dard deviation a for this model. 

A model which is more easily justified on physical grounds is the Weibull PDF, 

where S > 0, a > 0, and fi > 0. Despite the relevance of the Weibull distribution 14 to the 
strength of brittle materials, it is not often used, possibly because it is more difficult com- 
putationally to obtain reliability values than with the normal model. 
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If S follows the Weibull PDF, then ln(S) will have an extreme value distribution with PDF 
fmin (S) = ^-exp - exp (^1^] j ■ (5) 

The distribution of - ln(s) is 

. ( 6 ) 

Both of the above formulas are referred to as extreme value distributions. The use of ex- 
treme value distributions in a stress-strength model is illustrated in Figure 2. The extreme 
value distribution parameters are related to the Weibull parameters as follows: 

b = jj and u = - log a . 

In order to obtain the population Weibull shape and scale parameters /? and a from the 
known population mean ft and standard deviation a, the following approximations are suggested: 

p = 1.27 ft/< o - 0.56 (7) 

and 

a = P / T + 1 j . 

Stress-strength: extreme value 



3.0 3.S 4.0 

Figure 2. Stress-strength extreme value functions. 
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The functions defined in Equations 3, 4, 5, and 6 clearly have different shapes and they 
exhibit dramatically different tail behavior. Since reliability estimates depend strongly on the 
extreme upper tail of the stress PDF and the extreme lower tail of the strength PDF, the 
choice of model will typically have a substantial effect on the reliability estimate. For exam- 
ple, R is usually higher when calculated from the normal distribution than when the extreme 
value model is assumed. 

Applying PDFs Oat are capable of obtaining accurate high reliability estimates (e.g., 0.9(6)) 
requires prior knowledge of the functional form of the population PDF in addition to the 
availability of large data sets (e.g., 1,000 replicate specimens). For lower reliability values 
(e.g., 0.9), a goodness-of-fit test for PDF identification with a moderate amount of data is gen- 
erally adequate. The consequence of incorrect PDF selection and limited sample sizes are dis- 
cussed later in this report 


METHODS FOR COMPUTING RELIABILITY R 

In determining R from Equation 2 it should be noted that the integration process does 
not determine an area. The area A described by the intersecting functions in Figure 3 does 
not represent a 1 - R failure probability. The area A is the probability (P) that either S < T 
or s > T, that is, 

A « P(S < T) + P(s > T) , (8) 

where T is the point of intersection of the two functions. The area A is obviously not the 
same as the 1 - R from Equation 2 which determines P(S > s) jointly with P(s). 


StrMS-sfcangth ntfabitty modal: normal-normal 
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Numerical integration procedures are usually suggested if a closed form solution of 
Equation 2 is not available. The numerical integration process involves repeated application 
of a method such as Simpson’s Rule. The inner integral in Equation 2 is evaluated numeri- 
cally for each ordered S| value t * 1, ...; n resulting in an Ij(i) array of values. Each of 
Ij(i) is then multiplied by the corresponding f 2 (sj) forming another array I 2 (i) = f 2 (*i)Ii(i)- 
R is obtained from the I 2 array by reapplying the numerical integration method. This process 
will us ually provide accurate results for 51 ss n £ 101, where n is the number of mesh 
points in the integration process. Simulation results showed that the limits of integration can 
be obtained from ± six standard deviations from the mean. 

Closed form solutions are available when the assumed stress and strength PDFs are both 
normal or both Weibull. 

R Computation from Closed Form Solution 

If both stress and strength data can be represented by normal PDFs N(/t s , a t ) and 
N(/<s, as 2 )* respectively, then, 

R=P(S>s) "*®|)' (9) 

where (4>) is the standard N(0, 1) normal cumulative distribution function, n% and ps are 
means, and a s 2 and os 2 are variances of the stress and strength, respectively. 

If both fi and f 2 in Equation 2 are Weibull with different scale parameters ct\ and a 2 > 
but with a common shape parameter fi, then the integration indicated in Equation 2 gives the 
following closed form expression 12 



The common shape parameter means that both the stress and strength are skewed in the 
same way, which is a serious limitation. It is much more reasonable to have a stress distribu- 
tion with a heavy upper tail and a strength distribution with a heavy lower tail, but this is 
not possible unless the shape parameters can be varied separately. 

Nonpara metric Method 

THU method does not assume a PDF for either stress or strength data. It determines reli- 
ability from the ordered array of m stress (s) and n strength (S) values, where each of the S 
values are compared with all s values. R is the proportion of times S > s for the total num- 
ber of comparisons, that is 

R = Inn 2 2 a i* where a i = fl,Si>Sj (11) 

mn j-i i-1 jo, *j *S, 
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This method is not useful for obtaining high reliability even for relatively large data sets. 

It is obvious from Equation 11 that for high reliability calculations, mn must be very large; 
for example, 10* would be required in order to obtain R of 0.9(6). 

The Weibull, normal, or other parametric PDFs can provide estimates of high R values 
because of their ability to extrapolate beyond the available empirical data. Unfortunately, the 
amount of extrapolation dependency determines the magnitude of relative error in R. 

CONTAMINATED PROBABILITY DENSITY FUNCTIONS 

In order to illustrate the sensitivity of high reliability calculations to small deviations from 
assumed models, we will take the following approach. Consider the situation where with a 
high probability of 1 - e, specimens are obtained from a primary PDF, while with probability e, 
specimens come from a secondary PDF. Tim probability model is referred to as a contaminated 
model The secondary component is called the contamination, and the probability e is the 
amount of contamination. 

An example may help clarify this idea. Consider the situation where 97% of the time 
a specimen is obtained from a population of "good” specimens while the remaining 3% of 
the time consistently lower strength measurements are obtained, either due to manufactur- 
ing defects or to faulty testing. The primary PDF would correspond to the "good" speci- 
mens, the contamination would represent the distribution of flawed specimens, and the 
amount of contamination would be e = 0.03. 

The following procedure is introduced in order to examine the effects of computing 
high reliability values when uncertainties exist in selecting the functions for the stress- 
strength model. Initially, high reliability values are obtained from the normal stress- 
strength model (see Equation 9) using known PDFs with different mean values but equal 
coefficients of variation. The difference in mean values was determined from the required 
level of high reliability. Another R value is then obtained by applying this known distribu- 
tion with a small amount of contamination (e) in order to show an almost undetectable 
difference graphically between the true and contaminated PDFs. The effects of this differ- 
ence in the reliability computation are discussed in the following sections in order to exam- 
ine the sensitivity of the stress-strength model to the assumed PDFs. This procedure 
provides an effective way of demonstrating the effects of assuming a specific PDF in deter- 
mining high reliability. 

The normal PDF with variance contamination for the strength data is, 

Nsy <M& ^v 2 ) = (1 - e) N (><& as 2 ) + eN (j*s, Kios 2 ), (12) 

where and as 2 are the mean and variance for the uncontaminated normal strength distribu- 
tion, Kj is a scaling factor, and 100 £ is the percent contamination. 

The strength distribution with location contamination is 

N Si >d a s 2 ) = (1 - e) N a s 2 ) +£N^ S ±K 2 <7 S , as* 2 ) , (13) 
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where K 2 is a scaling factor for the mean /*s, and the sign determines which tail of the distri- 
bution is to be contaminated and as* 2 i* the variance on ps ± K^s- The location contami- 
nated PDF (see Equation 13) din provide reliability estimates to represent the potential of a 
secondary failure mode. Contamination of the stress distribution would be similar to that in 
Equations 12 and 13. It was not necessary to include contaminated distributions for both 
stress and strength in order to show substantial reduction in the high reliability estimates. 

The strength PDF contamination was sufficient. 

A linear relationship to obtain R for the reliability models when a combination of both 
contaminated and uncontaminated stress and strength normal PDFs can be written as, 

R = (1 -ci) (1 -c 2 ) Roo + ci (1 -e 2 ) Rio + *2 (1 -ei) R 01 + «i e 2 R 11 ( 14 ) 

where 100 £1 and 100 c 2 are percent contamination for the stress and strength distribution, 
and the Ry values are obtained for the case of variance contamination only; that is. 



and for location contamination, Rk^ would be 

(16) 


Equation 14 can be extended to include all conibinations of variance and location contamina- 
tion simultaneously, but it was not necessary for this sensitivity analysis. In Equation 16, if i, 
j = 0, then there is no contamination; for i, j = 1, then both stress and strength are contami- 
nated. For example, if there is contamination of variance of strength only, then 

R = (1 - c 2 ) R 00 + c 2 R 01 (I?) 

where 

R 00 — ® 


and 

R 01 = ® 


LOWER RELIABILITY BOUND 

A conservative lower bound on the reliability is introduced in order to protect against 
incorrectly identifying statistical functions in determining high R. The bound is obtained 
from a method proposed by Bolotin, 1 ^ and modified to employ the extreme value PDFs 


^Sq-^So \ 
°Si 2 + *So 2 


^Sq-^Sq ) 
-So 2 + *«o 2 ' 


- * (*•-««) 
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(see Equations 5 and 6). The method provides more conservative bounds than would be ob- 
tained from standard methods which are dependent on the assumed PDFs. The selection of 
the extreme value functions provides additional conservatism because of their heavier tails. 

The method is simple to use and is not restricted to any specified PDF. The reliability 
bounds are (see Figure 4), 

1 - W! W 2 > R > (1 - W!>(1 - W 2 ) (18) 

where [(1 - Wj)(l • W 2 )] represents the probability s < si and S > Si, which can be a 
somewhat conservative estimate. 


Stress-strength extreme value model: Bolotin R-bound 



Figures. Bolotin reliability bounds using extreme value functions. 

The lower bound is then, 

Rl > (1 - Wj) (1 - W 2 ) , (19) 

where 

W 2 =/" f 2 (s)ds and Wj * fj(S)dS 
for any choice of si = Sj . 
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GOOONESS-OF-F1T TEST 

The capability of determining desired PDFs from empirical data was investigated. The 
choice of PDF will be shown in the following sections to have a substantial effect on high re- 
liability computations, so it is important to examine model selection procedures. A statistical 
test 16 of goodness-of-fit was introduced in addition to graphical displays in order to select the 
desired PDFs. Empirical data used in the investigation was obtained by randomly selecting a 
relatively large number of values from a known normal PDF. A comparison of known contam- 
inated PDFs and the uncontaminated PDF is made with respect to the empirical values. 


RESULTS AND DISCUSSIONS 


Variance Contamination 

Shown in Figure 5a are reliability computation results and a graphical display of a normal/nor- 
mal stress-strength reliability model, where a 1% (e = 0.01) variance contamination was intro- 
duced and scaled by Ki = 4. The graphical display was obtained from application of 
Equations 12 and 13, where N&i s » <*s 2 ) » defined in Equation 3. The graph shows an almost 
undetectable difference between the contaminated and uncontaminated PDFs. This indicates 
that the choice of c and K are reasonable with respect to the potential differences between 
assumed and actual PDFs. However, the reliability values differ substantially (0%) versus 
0.998989). This implies that either one failure in a million or 1011 failures in a million is 
predicted depending on the selection of PDFs which can differ in probability values by less 
than 0.0005 in the extreme tail regions (see Figure 5b). Using "good* representative PDFs in 
the stress-strength model in predicting only a single failure will occur in one million opera- 
tions (e.g., number of flight hours) for R * 0.9(6) can result in a severe anticonservative esti- 
mate since for almost identical PDFs, 1011 failures per million could also be predicted. 

Normal stress-strength models 



Figure 5a. Stress strength normal functions with and without variance contamination. 
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Normal atreas-strength moctts with contamination 



Figure 5b. Detailed region d Intaraacting function* (sea Figures*). 

The accuracy of the high R estimate* depends on the level of precision in defining the ex- 
treme tail of PDFs. This requires selecting a PDF from a data set that accurately represents 
the known population function in the extreme tail regions with a probability difference of 
much less than 0.0005. Unfortunate^, this would require an unrealistically large data set. In 
current practice, if a very large data set is not available, then PDFs are selected from smaller 
sets with reliance on the functional representation in regions less than first ordered or greater 
than the largest value. 

The stress-strength procedure is quite effective for the range of R values between 0.5 
and 0.95 since usually in the extrapolation process, a small difference in the extreme tail prob- 
abilities values will not effect the required accuracy in R. Reliability results from uncontami- 
nated and variance contaminated (c — 0.05 and Kj = 5) PDFs showed no differences for a 
known R = 0.95. Unfortunately, in order to obtain high reliability, extrapolation into the ex- 
treme tail of the PDFs is required, thereby increasing the required level of precision neces- 
sary to distinguish between, for example, 0.998 and 0.9(g). 

In order to demonstrate the uncertainties in selecting specific PDFs from empirical data 
when computing high reliability values, the following displays are shown in Figure 6. In Fig- 
ure 6a, a plot is shown of the empirical normal cumulative density function (CDF) and the 
corresponding contaminated and uncontaminated normal distribution functions where the mean 
is 50 and standard deviation (SD) is 5, with sample size n = 100. Reliability values are also 
tabulated from the stress-strength model results using all six candidate functions. For exam- 
ple, R(3, 5) is the reliability obtained from variance contamination of 3% and a scale of 5 
for variance. A statistical goodness-of-fit test 16 that measures the relative differences in the 
tail region of the distributions was applied in addition to visual inspection in order to estab- 
lish if each function could represent the CDF of the ranked data. Results showed this to be 
true; see Figure 6b for the tabulated observed significance level (OSL) which shows in all 
cases OSL > 0.05, a requirement for the assumed function to be considered from the same 
population as the empirical data. 
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Figure 6a Goodness-of-flt: empirical versus functional normal CDFs. 
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Figure 6b. Lower tail of empiricai/normal distributions (see Figure 6a). 
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The results show that although each distribution fits the data quite well (see Figure 6 b), 
there is a large relative difference in R values: 0.9(6) for R(U.C) and 0.9957 for R(3, 5). In 
Figure 7, the results are similar to those in Figure 5. The variance contamination was 1% with 
a scale factor of 6 for both 05 and a t . Again, although the functions are similar, the relative 
reliabilities differ substantially (0.9(6) versus 0.9977197). As was the case in Figure 5, severe conse- 
quences could exist if R = 0.9(6) * assumed and the actual reliability was 0.9977197. This could re- 
sult in a number of premature failures, 2280 in one million, compared to the assumed one failure 
in a million. The results showed a low level of sensitivity to the selection of the factor Ki. 

Normal stress-strength models 



Stitts or StranQtfi 

Figure 7. ReliabUlty/normaJ functions with and without variance contamination. 

Location Contamination 

In Figure 8 , reliability computation results and a graphical display of the stress-strength mod- 
els are shown. The contaminated functions were obtained from 1% (e = 0.01) location contami- 
nation as defined in Equation 13 where K 2 = 4 and the (-) value is used for strength and (+) 
value for stress. The contamination in this case represents a secondary failure mode not consid- 
ered when assuming a specified function from the test results. For example, ignoring the possibil- 
ity that one in every 100 parts may have a tower strength level, say 4 standard deviations from 
the mean, can result in the reliabilities tabulated in the figure. That is, for the assumed correct 
model, R » 0.9(6), and the actual case where there was a tower strength level having one 
chance in 100 of occurring resulted in R - 0.999459. Figure 9 provides similar results to those 
in Figure 8 except there is a greater difference in reliability values 0.9(6) versus 0.991012 due to 
a greater shift (Kj = 6 ) in the mean value for the contaminated PDF. With a 1% contamina- 
tion this result is predictable since one in a hundred times a failure should occur because 
Ms - K 20 is less than the mean of stress value. The above figure shows die consequences of 
not being able to identify the correct function because of the inability to always detect a flawed 
component. The result is the determination of an overly optimistic reliability value when the 
true reliability could actually be orders of magnitude less. 



ProoaMiy 


Normal stress-strength models 
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Figure 8. ReliabHity/normal functions with and without location contamination. 


Normal stress-strength models 
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The results in Figure 10 are similar to those in Figure 5 except these were obtained from 
e = 0.03 and Ki = 3. If the estimated R = 0.9987350 is obtained from the empirical data 
and a higher R value is required (R - 0.9(g)), then a material with either greater strength or 
less contamination would be required. In order to obtain the required 0.9(g) from the origi- 
nal contaminated model, a mean strength of 87 is required (see Figure 11). The mean of 87 
requirement may not be acceptable to the designer, but this situation can occur if there is a 
substant ial amount of dispersion in the strength data resulting in a long-tailed PDF. The 
above situation shows when a potentially over-design situation could occur because of the in- 
ability to identify the correct PDF in the stress-strength model dim to inherent sensitivity and 
lack of information in the tail regions. This could prevent a good design from being ac- 
cepted if it is required that the assessment of the design be based upon reliability only. 


Normal stress-strength models 



Figure 10. Reiiabfllty/streaa-strength with and without contamination. 

The display in Figure 12 presents four possible reliability values for the case where the 
means and standard deviation are: stress (24, 2.4) and strength (51, 5.1). The result from 
the uncontaminated normal is 0.9(g). R = 0.995043 was obtained from the contaminated PDF 
application. Since, as was shown previously, the stress-strength model will often provide ei- 
ther relatively very high or low R values depending upon the chance selection of the PDFs. 

In order to compensate for the uncertainty in selecting the PDFIs for stress and strength data, 
extreme value distributions are introduced (see Figure 2) in the reliability computation. This 
resulted in R = 0.999045. Unfortunately, this did not provide a value lower than the contami- 
nated model result of 0.9950428. In order to obtain additional conservatism in the R esti- 
mate, a modification of a method by Bolotin is examined involving the determination of lower 
bound on R (see Figure 4 and Equation 19) in conjunction with the extreme value PDFs. 

The resultant lower bound estimate of 0.9796063 provides a significantly lower value than that 
of the contaminated model. This was also true for all contaminated models in this study. 
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Normal stress-strength models 



Figure 11. Increased strength requirements for high reliability. 
Extreme value stress-strength model 
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This lower bound estimate could provide some security in estimating R, although results 
may be excessively conservative for some practical applications. In Table 1, the distribution 
of R values as a function of the sample size is presented. R values were obtained from 
repeated application of the uncontaminated stress-strength model of Figure 5 using randomly 
selected, normally distributed samples. For a sample size of 10, R ranges from 0.9(6) to 
0.998417 indicating the instability associated with small samples. Higher order quantiles 
(e.g., 60%) were not included since they were all greater than 0.9(6). 


Tabtol. DISTRIBUTION OF R VERSUS SAMPLE SIZE 


Reliability R 

Sample Size 



10 

50 

100 

1000 

0.1 

0.996417 

0.900032 

0.999600 

0.999966 


(1583) 

(68) 

(20) 

(2) 

1 

0.999160 

0.996968 

0.999991 

0.999996 


(840) 

(32) 


(2) 

10 

0.999943 

0.999994 

0.999996 

0.999999 


(57) 

(8) 

(2) 

(1) 

25 

0.999994 

0.996996 

0.999699 

0.996999 


(6) 


(D 

(D 

50 

0.999999 

0.960969 

■ iT* ■ * r '■ 

0.996699 


(D 

(1) 

(D 

(1) 


( ) Corresponding Number of FaSures Per Million 


A sample size of 50 or 100 provides reasonable stability, and a sample of 1000 shows es- 
sentially no variability. The results from Table 1 show that for a sample of 1000, an estimate 
of R = 0.9(6) would be acceptable. This is not necessarily correct since results from the 
table only address the sample size issue which is independent of the uncertainties in the PDF 
selection process. There are two requirements for obtaining accurate high reliability values 
from the strength model: large samples (n > 1000) and knowledge of the population PDF. 
Reliability estimates of 0.95 are much less sensitive to the PDF assumption. If there is a sec- 
ondary failure mode due to occasional undetected poor manufacturing of the material or an 
unusually large load occurs that is not accounted for in the design process, then unknown 
lower reliability values (R < 0.95) can exist. 

CONCLUSIONS 

High reliability estimates from application of the statistical stress-strength model can vary sub- 
stantially even for almost undetectable differences in the assumed stress and strength PDFs. Spec- 
ifying high R values (e.g., 0.9(6)) for acceptable structural design can result in higher failure rates 
than anticipated if the assumed PDFs contain shorter tails than actually exist Over-design situa- 
tions can also occur when excessively long-tailed PDFs are applied to the stress-strength model. 

An effective method for identifying this nonrobust behavior involved application of contaminated 
and uncontaminated PDFs in the determination of reliability values. 
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A suggested method for obtaining a lower bound on the reliability estimate provided po- 
tentially overly conservative results but was effective in determing values that were lower than 
any of the R values computed for the contaminated models. 

The authors’ position regarding the computation of high reliability of 0.9(6) agrees with 
Breiman 1 ' who says "The probability of failure Pf = 1 x 10"® has an Alice in Wonderland 
flavor and should be banned from nonfiction literature.” It is therefore recommended that if 
high reliability calculations are absolutely essential, then the results should be subjected to a 
sensitivity analysis using contaminated distributions. High reliability values are meaningful only 
when these values are not substantially affected by an amount of contamination (e) consistent 
with the sample sizes, and a severity of contamination which is identified by engineering 
judgement. 
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